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Abstract 

We report the first observation of a striking history dependence of the struc- 
ture function of the vortex matter in the peak effect regime in a Nb single 
crystal by using small angle neutron scattering combined with in situ magnetic 
susceptibility measurements. Metastable phases of vortex matter, supercooled 
vortex liquid and superheated vortex solid, have been identified. We interpret 
our results as direct structural evidence for a first-order solid-liquid transition 
at the peak effect. 
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A current subject of wide interest concerns the existence of a vortex solid- 
liquid transition in type-II superconductors In addition to providing a 
model system for the study of melting and freezing, the vortex matter sys- 
tem offers unprecedented opportunities in studying the effects of quenched 
disorder on phase transitions. Recently, the magnetization jump in high- Tc 
superconductors YBa2Cu307_5 (YBCO), a widely accepted thermodynamic 
signature of a vortex solid-liquid transition and the peak effect, a well- 
known vortex-lattice softening anomaly in many low-T, li and high-T^ §^ 
superconductors in which the critical current exhibits a peak rather than de- 
creasing monotonically with increasing temperature, are found to occur at the 
same temperature [0,^]. However, there is no direct structural evidence indi- 
cating whether the underlying phase transition is solid-to-solid, solid-to-liquid, 
or even liquid-to-liquid in origin. Since quenched disorder is known to have 
important consequences for phase transitions |^Jl^, whether a solid-liquid 
transition can occur when random pinning is effective has broad implications 
in condensed matter physics. 

Small angle neutron scattering (SANS) pj|-15| is presently the most pow- 



erful technique for probing the vortex structure in bulk superconductors. For 
high-Tc superconductors, the magnetic penetration depth is typically quite 
long, yielding a very weak vortex signal in SANS experiments, which has so 
far prevented a direct structural study of the peak effect in high-Tc systems. 
Recent effort has therefore focused on low- Tc systems where one can detect 
the peak effect and measure the vortex structure using SANS in the same 
temperature and field regime |14,15|, and for these purposes Nb is experi- 



mentally the most favorable system because of the short penetration depth 
and the availability of large, ultrahigh quality single crystals. However, it 
has been controversial as to whether there is a vortex solid-liquid transition 
in Nb and the two previous SANS studies of this issue were not conclusive 



1 14, 15 1 . In the first study |lj], a characteristic temperature at which the az- 
imuthal width of the Bragg peaks starts to increase was interpreted as the 
vortex-lattice melting transition. In the second study |15|, a nearly isotropic 



ring of scattering was observed in the peak-effect regime, but was interpreted 
as a disordered solid rather than a liquid. In both studies, only field-cooled 
vortex states were measured. We report here striking hysteresis in neutron 
diffraction patterns of the vortex state in the peak-effect regime of Nb. The 
metastability of supercooled and superheated vortex states is directly demon- 
strated by applying a weak perturbation and observing the evolution of the 
SANS patterns. These data provide direct evidence for a first-order vortex 
solid-liquid transition. 

Our experiments were performed using the 30-m SANS instrument NG-7 
at the NIST Center for Neutron Research. The increased flux on the SANS 
instrument due to a new liquid-hydrogen cold source turns out to be impor- 
tant, but the key improvement provided by the present experiment is that 
the SANS and the ac magnetic response of the vortex array are measured si- 
multaneously in situ. Thus one can correlate the ordering in the vortex state 
directly with the macroscopic vortex dynamics. The coil for the magnetic 
ac susceptibility measurements can also be used to apply a small ac field to 
dynamically perturb the vortex system, allowing an in situ determination of 
the metastability of the vortex states. 

The sample is a Nb single crystal of 99.998% purity, a cylinder with a 
slightly uneven diameter from 1.316 cm at one end to 1.169 cm at the other, 
and 2.48 cm in length. The incident neutron beam has a mean wavelength 
A=6.0 A and a bandwidth AA/A=0.11. The neutron beam traverses through 
the central region of the sample, defined by a cadmium mask (0.7 cm in 
diameter), along the cylindrical axis which coincides with the three-fold sym- 
metric <111> crystallographic direction. The dc magnetic field is applied by 



a horizontal superconducting magnet along the same direction. The absolute 
accuracy of the measured sample temperature is it 0.20 K with a temperature 
stability better than ± 0.025 K. 

The peak-effect regime of the Nb sample is determined in situ by measur- 
ing the characteristic dip in the temperature dependence of the real-part of 
the ac magnetic susceptibility x'i^), as shown in Fig. 1(a) for H=3.75 kOe. 
The pronounced diamagnetic dip in x'(^) of the ac susceptibility corresponds 
to a strong peak effect in the critical current (or nonlinear conductance) in 
the sample fs],^. The onset, the peak, and the end of the peak effect are de- 
noted by To{H), Tp{H), and Tc2{H), respectively. (In YBCO, the end of the 
peak effect is still far below Tc2{H) (5|-p.) Fig. 1(b) shows the window of our 
experiment. The upper critical field i?c2(4.20K) = 4.23 kOe of this sample 



is higher than that of Ref. |]l^] but similar to that of Ref. |14|. The lower 
critical field line, separating the Abrikosov and Meissner phases, is estimated 
from the first appearance of vortex scattering in the SANS. 

For each (T, H) , we measure the vortex SANS patterns for different ther- 
mal paths. For low temperatures the vortex SANS images show sharp Bragg 
peaks with sixfold symmetry in agreement with previous studies ||Jl|,|ll, in- 



dependent of the thermal history. An example is shown in the inset of Fig. 1(b) 
for H=3.75 kOe and T=3.50 K. However, the vortex SANS pattern starts to 
show striking history dependence as the peak-effect regime is approached. For 
clarity, we define the field-cooled (FC) state as when the sample is cooled to 
the target temperature in a magnetic field, while the zero-field-cooled (ZFC) 
state is reached by cooling the sample in zero field to the target temperature 
and then increasing the magnetic field to the desired value. A field-cooled- 
warming (FC-W) state is when the system is cooled in field to a much lower 
temperature (~ 2 K) then warmed back to the final temperature. 

For the FC path, the vortex SANS patterns show nearly isotropic rings for 



Tp < T < Tc2 and broad Bragg spots for T < Tp. There is no clear sharpening 
in the Bragg peaks when Tq is crossed. Only at a lower T < Tp do the Bragg 
peaks become sharp. In contrast, for the ZFC and the FC-W paths, the sharp 
Bragg spots are observed for all temperatures up to Tc2- Shown in the top 
panel of Fig.2 are the ZFC and FC images at H=3.75 kOe and T=AAO K, 
which is just below ro(3.75kOe)=4.50 K. The images in the mid panel are for 
H=4m kOe and T=4A0 K, which is 0.10 K above rp(4.0kOe)=4.30 K. The 
intensities at the radial maximum for the mid panel SANS data are plotted 
in the lower panel. The sharp Bragg spots for the ZFC state indicate a vortex 
lattice with long-range-order (LRO) |jl^, while the broad spots for the FC 
state suggest a disordered phase with short-range-order. 

The orientational order of the vortex assembly can be quantified by the 
azimuthal widths of the Bragg peaks at the radial position of the intensity 
maximum. The azimuthal widths AO of the Bragg peaks can be obtained by 
fitting six Gaussian peaks to the data, for each (T, H) and path. Likewise, 
the translational order of the vortex assembly can be quantified by the radial 
widths AQ of the Bragg peaks. The Gaussian widths are plotted in Fig.3. 
Clearly, the ZFC states are more ordered than the FC states, translationally 
and orientationally, across the peak-effect regime. To determine the positional 
correlation of the vortex lines along the field direction, the rocking curves 
(Bragg peak intensity vs. the relative angle between the neutron beam and 
the applied magnetic field) of the ZFC and FC states at H=3.75 kOe and 
T=4.40 K are also measured. The rocking-curve width for the FC state is 
30% larger than that of the ZFC state, suggesting that the vortex lines are 
entangled [17| in the FC states, but nearly straight in the ZFC (and FC-W) 



states. 

Since a simultaneous broadening in radial and azimuthal widths is charac- 
teristic of a liquid (or glass) the hysteresis in both A6 and AQ suggests 



a first-order vortex solid-liquid (or glass) transition. A controversial issue is 
the location of the underlying equilibrium phase transition in the vortex mat- 
ter relative to the position of the peak effect PJl9|-21|. One interpretation 



|Q,^j2^ places the conjectured vortex solid- liquid transition at Tp, con- 
sistent with the recent experiments |^J8| in YBCO where the magnetization 
jump was found to coincide with the peak effect (the two are also related in 
other high- Tc systems ||2^). Another widely held view, is based on the clas- 
sical Lindemann criterion which would place Tm at Tc2{H) for Nb, provided 
the vortex-lattice elastic moduli are well-behaved (2^. In this scenario, the 
FC disordered phase seen here (as well as in |14,15|) is a supercooled liquid 



and the thermodynamic ground state is an ordered solid across the entire 
peak-effect regime. The third scenario places at or below the onset of the 



peak effect |21 



We find that it is possible to experimentally determine the ground state, 
and consequently the approximate position of Tm, of the vortex system. For 
this purpose, we use the ac susceptibility coil to apply a small ac magnetic 
field to shake the vortex assembly and use the SANS to observe how the 
structure of the system evolves. Even though the true ground state of the 
vortex system may not have been reached in the time scale of the shaking 
experiment due to random pinning, the evolution of the diffraction patterns 
leaves little doubt regarding the nature of the ground state. 

As shown in Fig. 4(a), after applying an ac field of 3.3 Oe at 100 Hz (1 
kHz gives similar results) for 10'^ sec, the ordered ZFC vortex lattice becomes 
completely disordered. Preliminary time-dependent data show that the Bragg 
peaks start to disappear within the first 10^ sec of the shaking experiment. In 
contrast, no measurable difference can be found for the disordered FC vortex 
states before and after the ac field is applied (for up to 10^ sec). By using 
the same approach, we find that the FC disordered states for T < Tp are 
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metastable and the ordered ZFC state is the ground state, opposite to that 
for T > Tp. In the T < Tp regime, the metastabihty is obviously stronger 
since a much larger ac field is needed to change the metastable state. We find 
that an ac field of 50 Oe (at ?a 0.1 Hz, using the superconducting magnet at 
non-persistent mode) can crystallize the disordered FC states at T < Tp. The 



shaking effects of an ac magnetic field were also observed in transport [24| 



and ac magnetization |25| in the peak-effect region in 2H-NbSe2, and were 
interpreted as a structural re-organization in the vortex matter, consistent 
with our direct SANS observations here. 

The metastable nature of the ZFC ordered state for T > Tp can also 
be observed in a field ramping experiment. For example, at a field ramp 
(increase) rate of less than 5 Oe/sec, the final vortex state (at T > Tp) is 
always ordered. In contrast, a ramp rate larger than 40 Oe/sec always results 
in a disordered state. However, the field ramping experiment alone cannot rule 
out a trivial possibility that the sample was heated to above Tc2{H) during 
the field ramping by the induced screening current. If this happens, after 
the field ramping stops, the sample cools back to the set temperature in field 
and the final state is actually a FC state. In the shaking experiment, the ac 
susceptibility of the sample is monitored and serves as an in situ thermometer 
to ensure that the sample temperature never fluctuates to above Tc2{H) during 
the entire SANS run (1 hour). 

Thus we conclude that, for T > Tp, the ordered ZFC vortex lattice is 
a superheated state and the ground state of the vortex system is a disor- 
dered vortex liquid, while for T < Tp, the ground state is a vortex crystal 
(or Bragg glass [0), and the disordered FC state is a supercooled vortex 
liquid. A thermodynamic phase transition must have taken place in the re- 
gion of the peak effect, T^ ~ Tp. This is consistent with the observations in 
YBCO where the magnetization jump was found to coincide with the peak 



effect [0,^]. In principle, one should also observe a discontinuous change of 



Qo when the vortex solid melts. Unfortunately, the expected shift in Qo 
due to the density change at the vortex solid-liquid transition is below our 
resolution limit, as we estimate here. From thermodynamic considerations 



H=Am kOe, Tm ^ Tp=4.30 K (rc2=4.50 K), dHm/dT= -0.846 kOe/K, and 
C66~ 0.6x10^ erg/cm'^ assuming a reasonable k ~1 and Lindemann number 
cl=0.1, the upper bound for the jump AB at melting is about 0.2 G. For 



a triangular lattice before melting Qo{B = H) = 2vrw-^,/^, the expected 



AQo ^ 3xl0"^(5o which is far too small to be resolved by the present SANS 
instruments. In Fig. 4(b), the radial peak position Qo remains unchanged after 
shaking-induced melting, consistent with the above estimate. 

Another interesting feature in Fig. 3 is that there appears to be a charac- 
teristic temperature below which the vortex state is insensitive to its thermal 
history. For example, at H=3.75 kOe, the FC, ZFC, and FC-W states at 3.50 
K are all ordered. The crossover temperature, t ~ 0.85 in Fig. 3, is in fact 
very close to the melting line identified previously (which was based on 
the measurements of the FC states). We interpret this temperature as the 
lower limit for supercooling. 

A closer examination of the diffuse scattering rings in Figs.l&2 suggests 
that they are not completely isotropic, indicating some short-range correla- 
tions. We attribute this short-range order to the coupling of the vortex liquid 
to the atomic crystal lattice, since experimentally one finds that in the Bragg 
solid phase, the specific orientation of the vortex lattice is always fixed to a 
certain crystallographic orientation of the sample. (This is the case for all 
low- Tc and high- Tc systems ||l3|.) In the supercooled state (top-right in Fig. 
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2), the highly correlated vortex liquid |14|, reminiscent of the hexatic liquid 
phase in 2D melting |26|, is also likely related to this coupling. 

We point out that our results imply the absence of superheating in con- 
ventional transport experiments with a large drive current. This solves a 



well-known puzzle [27,^ in which the history dependence of the nonlinear 
resistance always vanishes at Tp{H). Only with extremely low drive currents 



may one observe the subtle effects of superheating in transport [29|. The su- 
perheating of vortex solid contrasts sharply with that in regular bulk solids 



where surface melting prevents superheating |3C|. Surface melting of vortex 



solid is likely suppressed by the surface barrier |29| in type-II superconductors. 

In summary, we have observed a striking hysteresis in the SANS patterns 
of the vortex matter in the peak effect regime of a high quality Nb single 
crystal. The FC vortex states are disordered across the entire peak effect 
region and can be supercooled to below the onset of the peak effect. In 
contrast, the ZFC or the FC-W vortex states are ordered up to Hc2{T). The 
ordered ZFC/FC-W vortex lattices near Hc2{T) are identified as superheated 
solids which can be melted (disordered) by applying a small ac magnetic field. 
We interpret our SANS hysteresis data as direct evidence for a first-order 
vortex solid-liquid (glass) transition. 
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Figure Captions: 

Fig.l: Peak effect and the field-temperature {H — T) pliase diagram of 
Nb. (a) Representative ac magnetic susceptibility x'(^) vs. T for H(ic='i-75 
kOe (field-cooled). i?ac=3.3 Oe and /=1.0 kHz. HdcWHac- (For details on ac 
magnetic measurements, see [^,^.) Inset: the global H — T phase diagram 
for the Nb crystal used in this study, (b) Expanded view of the H — T phase 
diagram for Nb (the shaded box in (a)). Two representative SANS images of 
the field-cooled vortex states are shown (no background subtraction). 

Fig. 2: History-dependent neutron diffraction patterns. r=4.40 K for all 
images. The SANS images of the ZFC and FC vortex states for H=3.75 kOe 
(top panel: below the onset of the peak effect) and ff=4.00 kOe (mid panel: 
near the upper end of the peak-effect regime). The thick arrows indicate how 
the SANS images evolve after applying a small ac magnetic field (see text). 
The lower panel shows the intensity data at the radial maximum as a func- 
tion of the azimuthal angle 9 for the ZFC and FC SANS data (background 
subtracted) at i7=4.00 kOe and r=4.40 K. 

Fig. 3: Temperature dependence of the SANS widths, (a) The azimuthal 
widths of the Bragg peaks at Qo and (b) the radial widths vs. reduced tem- 
perature for ZFC, FC, and FC-W vortex states. The FC-W data points at 
3.75 kOe and 4.40 K are indicated by the crosses (and marked by arrows). 
The instrumental resolution in width is marked in (a). The lines are guides 
for the eye. The dashed lines indicate Tp{H). 

Fig. 4: The effect of an ac field, i/ac=3.3 Oe and /=100.0 Hz, on the 
ordered ZFC state for if =4.00 kOe and r=4.40 K (above Tp{H)). (a) I{Qo) 
vs. azimuthal angle; (b) the intensity I{Q) (averaged over the azimuthal 
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angle) as a function of the momentum transfer Q. The location of the radial 
maximum Qg is found by fitting I{Q) to a Gaussian. The lines are Gaussian 
fittings. 
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Fig.2 (color) Ling et al. 
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Fig. 3 Ling et al. 




_2 I I I I I I I I I I I I I 

60 120 180 240 300 360 

e (deg.) 




. -| I , \ , \ , \ , \ L \ , \ , \ , \ 1 

0.004 0.006 0.008 0.01 0.012 

Q (A-') 



Fig.4 Ling et al. 



